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1.  INTRODUCTION 


A  better  understanding  of  the  electromagnetic  properties  of  aqueous 
electrolytes  has  important  implications  for  oceanography  and  biology.  The 
purpose  of  this  study  is  to  consider  the  conductivity  relaxation  time  of  sea¬ 
water.  Potentially,  this  phenomenon  influences  the  electrical  properties  of 
seawater  above  5GHz.  Also,  the  dielectric  resonant  bands  in  the  far-infrared 
affect  the  electrical  properties  beyond  5  GHz,  and  their  effects  are  considered 
based  on  available  information  in  the  literature. 

Although  the  study  of  aqueous  electrolytic  solutions  has  a  long 
history,  much  yet  remains  to  be  understood.  Essentially,  the  electrical 
properties  of  solutions  of  high  concentration  [>0.1  Molarity  (M) ]  (i.e.,  sea¬ 
water)  are  not  completely  understood  and  are  typically  represented  by  empir¬ 
ical  models.  To  improve  the  present  understanding  of  seawater  conductivity, 
the  determination  of  its  relaxation  time  is  important  because  it  gives  informa¬ 
tion  about  the  mobility  of  the  charge  carriers. 

To  see  how  the  properties  of  mobility  influence  the  impedance  of 
electrolytes,  consider  the  standard  definition  for  conductivity  (Reference  1) 

N  N 

*  -  jEj  (e.in+)  -  »-.»-!>  *  ° 

a) 

where  n.  .  is  the  number  density  of  the  ith  +,-  charge  carrier,  e  is  the 

+  f  *  * t 

charge  of  the  ith  +,-  charge  carrier,  n  .is  the  carrier  mobility  of  the  ith 
+  ,-  charge  carrier,  and  N  is  the  number  of  solutes.  Since  most  electrolytes 
composing  the  ocean  are  1:1,  it  is  assumed  that 
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The  conductivity  for  a  single  1:1  solute  simplifies  to 

9  -  P(/<+  +  **_)  •  (2) 


A  simple  expression  for  the  carrier  mobility  can  be  obtained  from 
the  following  equations  (Reference  2) 


M  = 


<V> 

nr 


and 


m 


+  mg  <v>  =  eE 


(3) 


where  <v>  is  the  mean  velocity  of  the  charge  carrier,  m  is  the  mass,  g  is  the 
damping  constant,  e  is  the  charge  of  the  charge  carrier,  and  E  is  the  applied 
electric  field.  Let  <v>  and  E  be  time  harmonic  as  e^"t.  The  solution  of 
Equation  (3)  results  in  the  following 


M  = 


1  + 


j  - 
1  8 


Mo  = 


e 

mg 


(4) 


Substituting  this  result  into  Equation  (2),  the  conductivity  becomes 


<7  = 


o+ 


1  + 


j 

J  8+ 


*r,|F 


(5) 


It  is  assumed  that  the  damping  forces  will  be  similar  for  both  the  positive  and 
negative  charge  carriers.  Therefore, 

£+  =  £_  =  &  . 
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The  "g"  parameter  is  the  reciprocal  of  the  conductivity  relaxation  time.  This 
interpretation  is  easily  seen  from  Equation  (3).  When  the  applied  E-field  is 
turned  off,  the  charge  carriers  relax  back  to  equilibrium  and  the  mean  velo¬ 
city  becomes 


<v>  =  < v(t=0)  >  e~tg  . 

The  conductivity  relaxation  time,  1/g,  is  designated  "rc."  Thus, 


a  =  p(m_.  +  )  - - - 

0+  °“  1  +  j  “ 

1  g 


a  -r— 
O  1 


<0  T, 


where  is  the  conductivity  relaxation  time  related  to  mobility.  The  factor 
<rQ  represents  the  conductivity  according  to  the  Debye-Falkenhagen  theory 
(Reference  3)  and  the  remaining  factor  represents  the  correction  for  inertia. 

The  impedance  can  now  be  determined  for  a  parallel  plate  probe  as 
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where  i  is  the  length  of  the  sample  and  A  is  the  cross-sectional  area  of  the 
plate.  R  is  the  DC  resistance,  L  is  the  effective  inductance,  and  C  is  the 
capacitance  of  the  solution  between  the  plates.  Since  the  ions  cannot  move 
instantaneously  with  the  applied  field,  the  current  will  lag  the  voltage.  The 
conductivity  relaxation  time,  tq,  is  determined  by  computing  L/Rq. 

Impedance  measurements  of  this  type  should  be  used  instead  of 
microwave  cavity  measurements,  which  are  influenced  by  both  the  conductive 
and  dielectric  properties  of  the  measured  medium.  By  keeping  the  frequency 
of  the  applied  field  low  enough  (<  1  GHz)  so  that  the  real  part  of  the  di¬ 
electric  constant  is  constant,  the  imaginary  part  is  very  small  and  the  con¬ 
ductivity  high  enough  so  that  the  capacitive  (i.e.,  dielectric)  character  is 
shorted  out;  then  measurement  can  readity  be  made.  The  main  problem 
remaining  is  minimizing  the  stray  capacitance  and  residual  inductance. 

Section  2  describes  the  impedance  experiment  designed  to  measure 
the  conductivity  relaxation  time  of  seawater  and  related  electrocutes.  Section 
3  reviews  the  present  model  used  for  describing  the  electrical  properties  of 
seawater  and  discusses  the  effect  of  conductivity  relaxation  and  other  mecha¬ 
nisms  on  this  model. 

2  .  ELECTRICAL  IMPEDANCE  MEASUREMENTS 

The  major  thrust  of  this  project  is  to  measure  the  electrical  imped¬ 
ance  of  seawater-related  electrolytes.  A  special  conductivity  cell  was  de¬ 
signed  and  constructed  for  this  purpose.  The  following  describes  the  experi¬ 
mental  apparatus  used  and  presents  the  results  of  the  experiment. 


ivu'wwm 
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2.1  Experimental  Apparatus 

Figure  1  illustrates  the  experimental  apparatus.  The  impedance 
measurements  were  made  with  a  Hewlett-Packard  (HP)  multifrequency  LCR 
meter  model  number  4275 A,  which  operates  at  ten  different  frequencies  be¬ 
tween  10  kHz  and  10  MHz.  The  LCR  meter  was  computer  controlled  by  an  HP 
9845 A  computer.  A  conductivity  test  cell  was  connected  to  the  LCR  meter  in 
such  a  way  that  the  shields  of  the  4-terminal  network  were  connected  to  the 
shield  enclosing  the  test  cell.  As  shown  in  Figure  1,  the  current  passing 
through  the  cell  also  was  carried  by  the  shield  thus  canceling  the  magnetic 
field  generated  by  the  central  conductor.  This  reduced  the  residual  induc¬ 
tance  of  the  network. 

Figure  2  shows  the  unique  design  of  the  modified  Jones  cell  that 
was  used  as  the  test  cell.  The  cell  can  be  filled  and  drained  without  discon¬ 
necting  it  from  the  LCR  meter.  Also,  the  cell  can  be  temperature  controlled 
by  an  external  bath.  The  c  11  constant  is  17.2  (±0.15)cm  1.  The  tempera¬ 
ture  of  the  cell  is  monitored  by  placing  a  probe  inside  a  well  in  contact  with 
the  circulating  bath. 

The  accuracy  of  the  impedance  measurements  depends  on  the  LCR 
meter  and  background  effects  of  the  test  cell.  The  LCR  meter  does  open- 
and  short-circuit  measurements  of  the  test  network  and  measures  the  stray 
capacitance  and  residual  inductance.  The  values  are  stored  and  then  sub¬ 
tracted  from  the  actual  measurements.  The  open-circuit  measurement  is  done 
with  the  cell  empty.  The  short-circuit  measurement  is  done  with  the  cell 
filled  with  mercury.  In  this  way,  most  of  the  stray  capacitance  and  residual 
inductance  are  accounted  for.  Without  this  capability  these  measurements 
would  not  be  possible.  The  LCR  meter  reduces  the  stray  capacitance  and  the 
residual  inductance  to  some  nominal  level.  These  levels  are  measured  before 
every  experiment.  Figure  3  shows  the  accuracy  of  the  magnitude  and  phase 
of  the  impedance. 
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Figure  1  Experimental  Apparatus  (vector  meters  measure  magnitude  and  phase) 
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TEST  FREQUENCY 


Equations  In  table  represent 


Impedance  accuracy 
Phase  angle  accuracy 


a,  1/a:  See  Figure  A  Accuracy  Coeffi¬ 
cient  6raph. 

6  measurement  range: 

-180.000*  -  +180. 000* 

Display  counts  for  )Z|  and  e  (normal 


0  -  19999 


0  •  18000 


•  18000 


'Approximate  value  (unspecified). 


Number  of  significant  digits  displayed 
for  |Z|  and  8  depend  on  test  signal  lev¬ 
el.  range  and  frequency  (5  digits  max.). 

Accuracies  In  lined  areas  are 
unspecified. 


ACCURACY  COEFFICIENTS 


IIIMIIIIIUH 


400k*4MH2(C) 


lOfc'IOOk*  IMr  IOMHz(L*C ) 


20k-200k-2MHz(L-C) 


40k-400k-4MHZ(L ) 


•C-a-IZ|  display  counts  (M«*  resolution) 


L*C»1* |Z|  display  counts 


(  high  resolution  ) 


Figure  3  LCR  Meter  Accuracy 
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The  test  cell  changes  between  the  open  and  short  circuit  measure¬ 
ments  and  the  actual  measurement  of  electrolytic  solutions.  This  is  because 
of  the  double-layer  effect  at  each  electrode  and  the  difference  in  self¬ 
inductance  between  a  good  conductor  (mercury)  and  a  fair  conductor  (the 
electrolytic  solutions).  These  effects  cannot  be  nulled  out  by  the  LCR  meter 
since  they  are  not  present  during  the  open  and  short  circuit  tests.  The 
double  layer  effect  is  caused  by  the  electronic-to-ionic  conductor  interface. 
The  result  is  the  formation  of  a  leaky  capacitor  at  each  electrode,  as  shown 
in  Figure  4.  Also,  the  self-inductance  of  the  test  cell  is  much  larger  for  the 
electrolytic  solutions  than  for  mercury.  Calculations  for  an  infinitely  long 
conductor  indicate  the  difference  to  be  around  5  nH.  However,  end  effects 
will  certainly  increase  this  number. 

2.2  Results 

Measurements  on  carbon  composition  resistors  were  made  to  test  the 
system  background.  Since  the  resistors  are  electronic  conductors,  the  effects 
of  conductivity  relaxation  will  occur  at  frequencies  much  higher  than  ionic 
conductors.  Two  resistors  were  used,  one  with  a  resistance  of  228fi  and  the 
other  with  a  resistance  of  7570 fi.  This  range  of  values  of  resistance  is  close 
to  those  of  the  electrolytic  solutions  studied.  Tables  1  and  2  list  the  results 
of  the  carbon  composition  resistors.  The  format  is  typical  of  all  experimental 
results.  Background  measurements  are  made  after  the  open-  and  short- 
circuit  calibration  measurements  [a  zero  offset  adjustment  (ZOA)],  separately. 
Thus,  the  stray  capacitance  and  residual  inductance  of  the  network  are 
recorded.  The  stray  or  parasitic  background  capacitance  is  measured  at 
40  kHz  and  the  residual  or  parasitic  background  inductance  is  measured  at 
10  MHz.  The  measured  background  capacitance  and  inductance  as  listed  are 
typically  the  meter's  limit  of  sensitivity.  The  device  to  be  measured  is  then 
connected  to  the  test  network  and  approximate  values  for  the  inductance  and 
capacitance  are 
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obtained.  Then  the  computer  directs  the  LCR  meter  to  take  impedance  mea¬ 
surements  at  the  ten  available  frequencies  between  10  kHz  and  10  MHz.  In 
addition  to  the  frequency  and  the  magnitude  and  phase  of  the  impedance,  the 
computer  lists  the  approximate  applied  voltage  and  current,  and  a  corrected 
magnitude  and  phase  impedance.  The  measured  cell  capacitance  is  removed 
from  the  measured  impedance  to  give  the  corrected  impedance.  The  cell 
capacitance  in  the  case  of  the  carbon  composition  resistors  is  the  capacitance 
between  the  two  metal  contacts  at  opposite  ends  of  the  resistor.  In  the  case 
of  the  modified  Jones  cell,  it  is  typically  the  capacitance  between  the  two 
spherical  ends  of  the  capillary  tube.  For  the  carbon  composition  resistors, 
the  cell  capacitance  is  very  small  (0.4  to  0.1  pF  or  less).  The  impedance  of 
the  carbon  composition  resistors  is  largely  resistive  with  a  small  capacitive 
influence  at  the  higher  frequencies.  Notice  that  the  listed  approximate  induc¬ 
tance  is  negative.  This  means  the  impedance  was  actually  capacitive  and  not 
inductive.  Thus  the  LCR  meter  could  not  properly  interpret  the  result. 

Tables  3,  4,  and  5  list  the  main  results  of  the  experiment.  For 
sodium  chloride  and  seawater  solutions  (concentrations  are  listed  in  parts  per 
thousand  of  solute  over  solution  by  weight) ,  the  impedance  is  purely  resistive 
within  the  accuracy  of  the  LCR  meter.  The  inductive  effect  at  the  higher 
frequencies  is  less  than  0.1  percent.  In  Table  3,  the  cell  capacitance  domi¬ 
nates  at  all  frequencies  and  only  the  corrected  impedance  shows  an  inductive 
character  because  of  the  high  resistance.  This  result  is  different  from 
previous  less-precise  observations  that  showed  a  1-  to  2-percent  inductive 
effect  at  the  high  frequency  end  (Reference  4).  It  is  interesting  to  note  that 
the  electrodes  in  the  present  experiment  have  approximately  the  same  surface 
area  and  separation  as  the  electrodes  in  the  previous  experiment.  Because  of 
the  capillary  tube  in  the  modified  Jones  cell,  however,  the  cell  constant  is 
approximately  an  order  of  magnitude  larger.  If  the  observed  inductance  in 
the  previous  experiments  was  caused  by  some  electrode  process,  it  would 
explain  the  order  of  magnitude  smaller  inductive  effect  in  the  present 
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4  «  yE-04 

7.0OE-01 

3.00E-03 

2. 2995E+02 
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2. 2995H+02 

-0.02 

1  .  yE  +  05 

7. OAE-Ol 

3. 9OE-03 

2.2995S+02 

-0.02 

2. 2995E+02 

+  0.09 

2. OS +05 

7.00E-01 

3.09E-03 

2. 2992E+02 

-C.  Ol 

2 . 2992S  +  02 

+  0. 03 

4 ,  C  E  *  o  5 
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+  0. 00 
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3. 0OE-O3 

2.2953E+02 

♦  0.05 

2.2953E+02 

+  0.39 

4 .  C  E  ►  S': 
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2. 2918E+02 
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+  0. 74 

1 . 0E+07 
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2. 293OE+02 
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+  1 .  U  9 
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++ 

+  +  +  ■♦•  +  +4-  f  «f  • 

C++++++++H 

MEASURED  IMPEDANCE 

•r4+  +  +  +  Tf+  +  4+>+!‘t  +  +  + 

CORRECTED  IMPEDE. UC 

r  R  E  G  ;J  E  ?l  C  Y 

’  VOL “AG 7 

CUPPEMT 

HHGMITUI-E 

PHASE 

MAGN I  Tuts 

P  h  .9  3  r* 

:  H  2 ) 

(VOLTS) 

(AMPS) 

<  OHMS) 

v  I'  E  u  ■' 

'  Oris  ,• 

I'  c  C  > 

1 .  OE  '•  •  1 

7 . OOE-O 1 

3. 00E-O3 

2. 2957 -02 

— 0 .  1  1 

2. 2957E+02 

-o.  10 

2. 0E  +  0- 

7. COE-01 

3. COE-03 

2.29512*02 

-0.05 

2. 2951E  +  02 

-0.  ,  5 

4 . OR  +  04 

7. 0OE-01 

3. O0E-O3 

2. 2947E+02 

-0. 03 

2.2+47E-02 

-0.02 

1 . 0E  *  95 

7. OOE-C1 

3 . 0  0  E  -  ©  3 

2 . 2945E+02 

-0.02 

2.29455*02 

-  0  .  C*  0 

2. 0E+05 

7.O0E-01 

i . O0E-03 

2 . 2942E+02 

-0.01 

2  ■  i’M  i'E  *  C  — 

♦  0 , 0  9 

4. 0E+ v? 

<2 . 9OE-0 1 

3. 0OE-O3 

2 . 2937E+02 

+0.01 

<l  *  ©  2 

+  0 . 0  7 

1 . 0  z  * »' 

*£ .  ?  G  £  *•  0 1 

3. 90E-C3 

2.2933E+02 

+  0. 00 

2 . 2932E  ►  0 2 

*0 .  1  £ 

Xi  •  02 J  0  £ 

7. 00E-01 

3.  OC'E-OS 

2. 2923E+02 

+  o .  o  ? 

2 . 2  9  2  3  E  +  0  2 

K‘ .  37 

4 ,  OE  ‘  0 o 

i. 9OE-01 

3. COE-03 

2. 2S94E+02 

♦0.07 

2. 2892E+02 

+  0. 7  1 

1 . Cc+07 
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Table  5 

SIMULATED  SEAWATER  SOLUTION  IMPEDANCE 
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FREOUSNCY  PARASITIC  BACKGROUND  CAPACITANCE  AND  CONDUCTANCE 
4C!  6 . 2COOO0OO0OOE- 1 ?  -.000000000023 


FREQUENCY  PARASITIC  BACKGROUND  INDUCTANCE 
1  011  6 . 325O0000000E-  1 0 


T=23.SC 


ARRkOXI.IATE  INDUCTANCE 
,00  0 .0  0 0 ,0  1294 


CELL  CAPACITANCE 
1 . 24000000000E- 12 


MEASURED  IMPEDANCE  CORRECTED  IMPEDANCE 


FREQUENCY 

VOLTAGE 

CURRENT 

MAGNITUDE 

PHASE 

MAGNITUDE 

PHASE 

'•HI) 

< VOLTS) 

<  AMPS ) 

<  OHMS  > 

(DEG) 

(OHMS) 

(DEG) 

1 . OE  +  04 

7.30E-01 

2.20E--03 

3. 4419E+02 

-0.  09 

3.4419E+02 

-0.09 

2. 0E  +  04 

7.30E-O1 

2.20E-33 

3.4417E+02 

-O.05 

3 . 44 1 7E+02 

-0.04 

4.0E+04 

7.30E-01 

2.20E-O3 

3. 4415E+02 

-0.  03 

3. 4 4 1 5 e  +  0 2 

-0.O2 

1 . OE+05 

7 . 30E-0 1 

2.20E-03 
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-0.O2 
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-0.00 

2. 0E  +  03 
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4. 0E  +  C5 

7. 7OE-01 

2.2OE-03 
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+  0.01 

3. 4407E+02 
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1 , 0  E  + 0  5 

7.7OE-01 

2. 20E-03 

3 . 4403E+O2 

-0.01 

3 . 4404E+C2 

+  0.13 

2. 0E+05 

7. 70E-O1 

2 • 20E— 03 

3.4392E+02 

+  0. 0'» 

3.4391E+02 

+  0 .  34 

4  >  0E+&S 

7.  50E-01 

2.20E-03 

3.4351E+02 

+  0 . 04 

3. 4343E+02 

+0. 5  3 

i . oe+ot 

6. 7OE-01 

1.90E-33 

3.4398E+02 

+0.13 

3.  4-383E+02 

+  1.37 

4.  a?  r 
♦  ♦  •**  * 

ftPP * OX  I  MATE  I N&UCTANCE 
•  C*  *5  U  ,0#v  1  *,5  *4 


CEl.L  CAPACITANCE 
1 .  jOOOOv.  >OO0OE-  12 


♦  +  ♦  *r  ♦  H 


MEASURED  IMPEDANCE  CORRECTED  IMPE2AHCE 


r  PSOUE'.CY 

VOLTAGE 

CURRENT 

MAGNITUDE 

PHASE 

MAGNITUDE 

PHASE 

1  nZ  1 

•VOLTS) 

(AMPS) 

(OHMS ) 

'DEG) 

(OHMS > 

<  I.  E  G  > 

1 . 0E+04 

7.3OE-01 

2.2OE-03 

3.441 3E+02 

-0.  09 

3. 4413E+02 

-0. 09 

2. 0E+04 

7. 3OE-01 

2. 20c -03 

3.441 12+02 

-0 . 04 

3.441 1E+02 

-  0 . 0  4 

4. OS-04 
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3 . 4397E  >02 

♦  0 . 0  3 

4.  VE  -C  3 

7. 70E-01 

2.20E-03 

3. 4339E+02 

+  0 . 0  1 

3. 433? E +02 

o .  o  r 

1.CE+0S 

7. 70E-0 1 

2 . 20E-03 

3. 43342+02 

-0. 01 

3. 43042- 0 2 

-0 .  1  5 

»  -  «  0  £  ♦  -  “* 

7. 702-0 1 

2.20E-03 

3. 4359E+02 

+  0 . 04 

3. 4  3532  +  02 

-  c* *  s  £ 

4.0= - O 5 

7. 502-31 

2.20E-03 

3. 4327E+02 

+  0. 04 

3. 4  324E  +  02 

*  0  .  £  <?• 

l . 02  +  07 

5. 70E-01 

t . A0E-O3 

3. 4372E+02 

+  0.13 

3 . 4  3  3 5  E  ♦  0  2 

♦  1 . 74 

I* 

J4 
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experiment.  That  is,  the  electrode  process  remained  unchanged  while  the  cell 
impedance  was  increased  by  an  order  of  magnitude.  Further,  the  relative 
changes  in  the  impedance  are  essentially  the  same  in  the  previous  and  present 
experiments . 

Another  set  of  measurements  attempted  to  observe  the  Debye- 
Falkenhagen  effect  (Reference  3)  that  is  observable  within  the  frequency 
band  of  the  LCR  meter  but  at  ve^  low  concentrations  (around  1000  times  less 
than  seawater  concentrations).  Tables  6  and  7  show  the  results  for  0.01- 
Molarity  (M)  solutions.  The  Debye-Falkenhagen  effect  is  much  stronger  for 
MgSC>4  than  for  KC1;  thus,  a  comparison  can  be  made  between  the  two  tables. 
Within  the  accuracy  of  the  measurements,  no  differences  are  observed  and  the 
corrected  impedance  is  purely  resistive.  The  concentration  is  too  high  to 
observe  this  effect. 

Finally,  much  of  the  initial  motivation  to  carry  out  these 

experiments  was  to  measure  the  conductivity  relaxation  time.  This  number  is 
-13 

at  least  10  sec,  which  is  based  on  the  initial  assumptions  of  Debye  and 
Falkenhagen  (Reference  3).  Based  on  the  experimental  results  of  Table  5, 
the  maximum  seawater  conductivity  relaxation  time,  r  ,  is  5  x  10  sec. 
Therefore , 

1  x  10  13  <  rc<5  x  10  1(*sec  or 
3.2  x  108  (  f  (  =  x— i — )<  1.6  x  1 0 1 2Hz 

N  C  2  ir  r  c 

3.  PRESENT  UNDERSTANDING  OF  THE  ELECTRICAL  PROPERTIES  OF 

SEAWATER 

The  problem  of  understanding  the  electrical  properties  of  seawater 
is  a  very  complex  one.  This  is  because  the  detailed  physics  of  electrolytic 
solutions  is  not  completely  understood  and  seawater  is  a  heterogeneous  oolu- 
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tion  of  many  electrolytes.  The  importance  of  the  electrical  properties  (i.e., 
dielectric  and  conductive)  at  radar  frequencies  has  required  some  form  of 
tabulation  as  a  function  of  temperature,  salinity,  and  frequency.  This  is 
currently  done  in  a  semi-empirical  way  using  the  formulas  of  Debye,  and  Cole 
and  Cole  (References  5  and  6)  fitted  to  experimental  data  (Reference  7). 
The  result  is  a  reasonable  representation  of  the  true  dielectric  and  conductive 
values  of  seawater  as  a  function  of  temperature  and  salinity  below  5  GHz. 
However,  precision  microwave  radiometry  requires  equal  precision  of  the 
known  electrical  properties  of  seawater.  To  improve  the  present  semi- 
empirical  formulas  requires  a  greater  understanding  of  the  theory,  including 
observed  anomalous  electrical  properties.  A  brief  review  of  the  present 
understanding  of  the  electrical  properties  of  seawater  and  seawater-related 
electrolytes  is  presented  in  the  following  section. 

3.1  Microwave  Electrical  Properties 


Microwave  data  of  the  complex  relative  permittivity  for  electrolytic 
solutions  is  typically  fitted  to  an  equation  of  the  form  (Reference  7) 


( u>  •  T ,  S )  -  e„,  * 


's(T,S)  -  «, 

1  +  (j«rn(T,S)? 


_  i  g-  (T , S )  /o\ 

-  j  - - (8) 


where 


temperature 
sal  ini ty 

the  relative  dielectric  constant  at  »»  l/r^ 
the  static  relative  dielectric  constant 
permittivity  of  free  space  [Farads/mJ 
2irf,  f  ■  frequency  in  hertz 
conductivity  [S/m] 

Debye  relaxation  time 

empirical  spreading  parameter  of  relaxation  times 
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The  conductivity  relaxation  process  is  not  represented  by  Equation  (8)  and 
neither  are  higher  frequency  resonances  which  impact  .  However,  since 
this  model  only  claims  agreement  below  5  GHz,  the  details  of  the  Debye  and 
conductivity  relaxations  are  not  required.  For  example,  for  «  =  0  Equation 
(3)  can  be  expanded  to  the  form 

eU,T,S)  =  +  (*g  -  eM)  (1  -  (i.rD)  -  (,»rD)2-  — )  -  j  -J-  (9) 

-  *,  -  <«8  -  >»  )  <-D>2  -  it(«s  -  ‘„><-’D>  *  ^1 

Because  «  1  below  5  GHz,  this  expansion  will  work  very  well.  The 

expansion  of  any  proper  function  would  produce  the  sane  first-order  struc¬ 
ture.  Thus,  curve  fitting  data  to  Equation  (8)  at  frequencies  below  5  GHz 
will  not  definitively  determine  the  true  structure  of  the  dielectric  properties 
beyond  5  GHz.  This  fact  is  recognised  by  Klein  and  Swift  (Reference  7). 

In  principle,  and  <r  are  functions  of  temperature,  salinity,  and 

frequency,  and  and  rD  are  functions  of  temperature  and  salinity.  A  more 

general  form  of  Equation  (9)  is  given  by 


*  ( ui ,  T ,  S )  *  *  o  ( u> ,  T ,  S )  + 


*  (T,  S)  -  *,„(«, T,S) 


1  +  [ r  d(T,S) ] 


TT75T 


(10) 
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where 


<tD£  =  DC  conductivity 
r  =  conductivity  relaxation  time 


The  Debye  relaxation  process  accounts  for  the  damped  rotation  of  the  water 
molecule  in  the  liquid  state.  The  conductivity  relaxation  process  accounts  for 
the  damped  translational  motion  of  the  ionic  charge  carrier.  Debye  in  1913 
(Reference  8)  approximately  determined  the  Debye  relaxation  time,  by 

considering  Stoke's  law  converted  for  a  rotating  sphere  and  Brownian  motion. 
The  following  formula  was  obtained. 


r 


D 


birq  8^ 

“ET 


(11) 


where 

v  -  viscosity  of  water 
a  =  radius  of  water  molecule 
k  =  Boltzmann's  constant 


Surprisingly,  this  formula  predicts  reasonably  well  the  observed  relaxation 
time  of  water,  which  is  on  the  order  of  1  x  10  11  sec.  Intuitively,  the 
relaxation  time  of  the  translational  motion  of  the  ionic  charge  carriers  in  water 
should  be  similar  to  the  rotational  relaxation  of  the  water  molecule.  Following 
a  similar  approach  as  Debye,  the  conductivity  relaxation  time  can  be  approxi¬ 
mately  found  by  examining  the  charge  density  in  light  of  Stoke's  law  and 
Brownian  motion.  The  result  is 

6wri  bx“  (12) 
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where 


b  =  radius  of  charge  carrier 


x  =  1/e  point  of  Maxwe 1 1 -Bo i t zmann  distribution 
o 


of  charge  density 


Since 


b  ~  a 


xo  >  a 


>  rD 


This  result  crudely  verifies  our  intuition.  It  is  interesting  to  note  that  the 
Debye-Falkenhagen  computation  of  the  conductivity  relaxation  time  is  based  on 
the  velocity  or  mobility  part  of  the  conductivity  and  does  not  include  the 
random  motion  of  the  ion. 


The  real  and  imaginary  parts  of  the  general  complex  permittivity 
[Equation  (10)]  are  for  «  =  0,  thus 


(13) 


e  =  + 
00 


e  _  -  e  „ 

S  00 


"DC  Cleo 


1  +  (o»Tjj)^  1  +  ( <0  r 


(14) 


e~-  (ts  ~  e<x>  )“TP 


"DC  TC/%W 


(15) 


1  +  (lUTjj) 


1  +  (u,TC)" 


■  -  V'.\v  ,'.y-  .* ■  »" ' 


•Lviviv: 


STD-R-1071 


The  Jotim  Hophtm  UnHranity 
APPUIO  PHYSICS  LAPOPATOPY 


L aural,  Maryland 


Figures  5a  and  5b  are  plots  of  «'  and  * " ,  respectively,  for  different  values 

of  t  .  The  cases  for  rc  =  or  0  are  the  same  curves.  Interestingly,  with 

t  =  td  the  conductivity  effects  cannot  be  differentiated  from  the  case  of 

r  =0.  However,  for  r  /t_  =  0.5  and  2  there  are  observable  differences.  If 
c  CD 

rc  ~  td  (it  is  unlikely  that  they  are  exactly  equal)  then  the  conductivity 
relaxation  effect  is  an  important  term  in  the  complex  permittivity  formula. 
The  efforts  to  model  the  permittivity  beyond  5  GHz  would  then  be  helped  by 
including  this  effect. 

Figures  6a  and  6b  illustrate  the  dielectric  properties  of  pure  water. 

Figure  6a  also  shows  typical  skin  depth  attenuation  of  seawater  (References  9 

and  10).  The  pure  water  spectra  show  a  fair  amount  of  structure.  The 

hump  at  20  GHz  is  the  Debye  relaxation  phenomenon.  The  two  peaks  between 
12  13 

10  and  3  x  10  Hz  are  called  the  hindered  translation,  v^,  and  librational 
band,  i^,  respectively.  They  are  caused  by  H90  and  HgO  (i.e.,  dimmer¬ 
like)  interactions  (Reference  11).  Most  of  the  higher  frequency  peaks  are 
caused  by  the  internal  motion  of  the  water  molecule  (i.e.,  rotation  vibration 
(Reference  11).  Beyond  the  Debye  relaxation,  the  complex  permittivity  is 
represented  by  .  It  is  dominated  by  these  infrared  resonances.  This 
can  be  seen  by  examining  the  real  part  of  £,*,(=  j  ); 

4.9  >  t'n  >  1.69  >  1 
Microwave  Optical 

is  determined  by  resonances  and  relaxations  beyond  the  frequency  of 
interest.  Based  on  the  above  relationship,  the  microwave  region  is  heavily 
influenced  by  the  infrared  resonances. 
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The  presence  of  a  solute  does  influence  «,*, .  This  is  demonstrated 
by  Figure  7,  which  shows  experimental  data  of  the  and  v,j,  band  of  pure 
D20  and  DjO  KC1  and  NaCl  solutions  at  4  Molarity  (M).  The  solutions  are 
highly  concentrated  compared  to  seawater  but  do  show  significant  changes  in 
the  spectrum .  This  will  directly  impact  at  microwave  frequencies  and  make 
it  a  function  of  salinity.  It  is  also  well  know  (Reference  11)  that  these  bands 
are  sensitive  to  temperature.  Thus  is  indeed  a  function  of  frequency, 
salinity,  and  temperature. 

The  relaxation  time  spreading  parameter,  a  ,  is  not  zero  and  is 
thought  to  be  a  function  of  temperature  and  salinity  (References  12  and  13). 
This  is  illustrated  in  Figure  8  for  HgO  solutions  of  NaCl.  The  precise  nature 
of  a  will  not  be  understood  until  the  conductivity  relaxation  effects  are  un¬ 
derstood.  The  relaxation  time  spreading  parameter  influences  the  falloff  of 
the  Debye  relaxation.  As  illustrated  by  Figures  5a  and  5b,  the  conductivity 
relaxation  process  also  influences  the  falloff  of  the  curves.  Thus,  these  two 
different  yet  similar  effects  may  be  difficult  to  separate. 


3.2  Anomalous  Electrical  Properties 

In  addition  to  well-understood  physical  mechanisms  that  impact 
seawater  electrical  properties,  there  are  reported  observations  of  anomalous 
behavior  of  solutions  of  single  electrolytics  which  compose  seawater.  These 
are  listed  below  (Reference  13). 

(1)  Anomalously  high  conductivity  of  electrolytic  solutions  at 
concentrations  close  to  seawater  at  582  MHz  (Reference  14). 

(2)  Additional  relaxation  phenomena  at  submillimeter  wavelengths. 
The  dielectric  constant  is  modelled  as 


4-Molarity  Solutions 


Figure  7 


Experimental  Data  of  the  and  vT  Band 
DgO  and  DjO  KC1  and  NaCl  Solutions  at  4 
Concentration 


of  Pure 
Molarity 
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1  +  (  j«>  t) 


where  is  the  additional  relaxation  time. 


4.  CONCLUSIONS 

The  electrical  pro"  arties  of  seawater  have  been  examined  with 
particular  emphasis  on  the  effects  of  relaxation  of  the  charge  carriers  with 
the  frequency  of  the  applied  field.  Such  effects  have  not  been  definitively 
quantified  either  experimentally  or  theoretically.  Experiments  conducted  were 
not  conclusive.  It  is  possible  that  the  conductivity  relaxation  time  (the 
characteristic  time  its  takes  a  translating  charge  carrier  to  relax  back  to 
equilibrium)  is  similar  to  the  Debye  relaxation  tine  (the  characteristic  time  it 
takes  a  rotating  molecule  to  relax  back  to  equilibrium  in  a  fluid  or  solid), 
making  direct  observation  of  this  effect  difficult.  The  conductivity  relaxation 
time  has  been  constrained  to  be  within  5  x  10  ^  ^  >  1  x  10  ^  sec.  It  is 

quite  possible  that  ~  Tq  since  the  dampening  mechanisms  should  be  similar. 
If  this  is  true,  the  dielectric  models  will  need  to  include  the  effects  of 
conductivitj’-  relaxation.  Also,  it  is  recognized  that  resonant  far-infrared 
(including  submillimeter)  bands  can  influence  the  dielectric  properties  at  and 
above  40  GHz. 

A  more  definitive  impedance  experiment  attempting  to  measure  the 
conductivity  relaxation  time  will  require  the  following: 

•  A  test  cell  with  a  much  larger  cell  constant 

•  Impedance  measurements  at  frequencies  close  to  1  GHz 

•  A  high  degree  of  precision  or  a  means  of  nulling  the  DC 
resistance  of  the  electrolytic  solution 

Also,  the  development  of  a  reasonable  theory  may  aid  in  the  selection  of  a 
more  optimal  temperature  and  salinity  of  the  samples. 
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